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Abstract 

We calculated the electron mobility of 14 two dimensional semiconductors with composition 
of MX 2 , where M (= Mo, W, Sn, Hf, Zr and Pt) is the transition metal, and X is S, Se and 
Te. We treated the scattering matrix by deformation potential approximation. Long wave lon¬ 
gitudinal acoustical and optical phonon scatterings are included. Piezoelectric scattering in the 
compounds without inversion symmetry is also taken into account. We found that out of the 14 
compounds, WSe 2 , PtS 2 and PtSe 2 , are promising regarding to the possible high electron mobil¬ 
ity and finite band gap. The phonon limited mobility in PtSe 2 reaches about 3000 cm 2 V -1 s -1 
at room temperature which is the highest among the compounds. The bandgap under the local 
density approximation is 1.25 eV. Our results can be a guide for experiments to search for better 
two-dimensional materials for future semiconductor devices. 
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I. INTRODUCTION 


Two-dimensional (2D) layered materials have received a lot of research activities since 
the discovering of graphene. New materials such as transitional metal dichalcogenides, espe¬ 
cially MoS 2 [I], an d black phosphorus |2] with one or several atomic layers were synthesized. 
Field effect transistors, photo detectors and light emitting diodes were fabricated to demon¬ 
strate the potential applications of these materials. High speed radio frequency devices are 
fabricated which make full use of the ultrahigh electron mobility in graphene [3]. Carrier 
mobility is one of the key parameters for semiconductors to be used in the high speed or 
high frequency devices. Graphene has ultrahigh mobility, but it is intrinsically metallic. 
It is possible to open a gap at the Dirac cone, but usually the gap is tiny, and required 
extra parameters, for example, external electric fields^], which is not easy to fully integrate 
it with present semiconductor processes. When it was used in logical devices, the on/off 
ratio is small, and the power consumption is large at the off-state. Monolayer M 0 S 2 has a 
direct band gap about 1.7 eV. However, the electron mobility at room temperature is only 
up to ten cm 2 V -1 s -1 . There are quite a lot of efforts to improve the carrier mobility. The 
mobility can be enhanced by even 20-fold to ~200 cm 2 V _1 s _1 when encapsulated with a 
liigh-fc dielectric layer [5], This mobility enhancement was attributed to screening of the 
charged impurity scattering by the dielectric layer[2] • On the device level, the mobility can 
be engineered by electron doping. WSe 2 monolayer with a direct band gap about 1.65 eV 
was explored after M 0 S 2 IIJ, which shows high field-effect mobility about 140 cm 2 V _1 s _1 . It 
is already substantially higher than that of most of the reported room-temperature mobility 
values of M 0 S 2 . lH-MoTe 2 has a direct band gap a little above 1.0 eV. The room temper¬ 
ature electron mobility was found to be about 30 cm 2 V _1 s _1 without gatingP, 9j. This is 
already comparable with that of M 0 S 2 . 

According to the work of Yoon [TO]. due to the heavier electron effective mass and a 
lower mobility, M 0 S 2 is not ideal for high-performance semiconductor device applications. 
A natural question arises: can we find 2D materials with larger carrier mobility? And how 
far can we go to improve the mobility in the reported materials? Both of these questions 
can be answered by theoretical calculations. The ultramost mobility can be reached is lim¬ 
ited by the intrinsic scattering. Phonon scattering is the one which cannot be avoided at 
temperatures other than zero. There are quite some works to calculate the phonon-limited 
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carrier mobility in 2D materials. The value in M 0 S 2 was first calculated by Kaasbjerg et 
al CH The contributions from acoustic and optical phonons are included and electron- 
phonon coupling matrices are calculated by the frozen phonon method. Electron-phonon, as 
well as piezoelectric interactions, are taken into account. The calculated room temperature 
mobility is about 410 cm 2 V -1 s -1 which sets the upper bound of electron mobility. This 
value estimated by Restrepo et al. [13] is 225 cm 2 V" x s T Li[T4] et al. developed the method 
to calculate the electron-phonon interactions based on the density functional perturbation 
theory and obtained similar results. The phonon limited electron mobility in single-layer 
M 0 S 2 is dorminated by the acoustic phonons at low field, while the degeneration by the 
optical phonons happens at high field according to the results of Monte Carlo simulations 
using the corresponding deformation potentials [15]. Very recently, Kim et al. [16] calculated 
the phonon limited electron and hole mobility of four monolayer transition metal dichalco- 
genides: MoS 2 , MoSe 2 , WS 2 and WSe 2 . WS 2 provides the best performances of electrons 
and holes with high mobilities. It was also found that p-type WSe 2 shows hole mobilities 
comparable or even larger than that of bulk silicon at room temperature. 


Our previous work shows that MoTe 2 , ZrSe 2 and HfSe 2 have the largest electron mobility 
if only long wave acoustical phonon scattering was taken into account im However, as 
shown in the work of Kaasbjerg [IT], the optical phonon may also play equally important 
role in M 0 S 2 in the room-temperature range due to the low optical phonon frequencies. In 
this work, we take these effects into account, and try to screen out MX 2 type 2D monolayer 
semiconductors with possible high electron mobility. We performed electronic calculations 
of the selected semiconductors with composition of MX 2 , where M (= Mo W, Sn, Hf, Zr or 
Pt) is the transitional metal, and X is S, Se or Te. I 11 order to fast screen out the materials 
with high performances, mobility limited by long wave phonons was estimated by calculating 
the corresponding deformation potentials. We found that out of the 14 compounds, WSe 2 
with lH-structure and PtSe 2 with lT-structure, are promising regarding to their mobility 
and sizable band gap. We can see that due to the optical phonon scatter, only PtS 2 and 
PtSe 2 with the lT-structure can have electron mobility larger than 1500 cm 2 V _1 s _1 . 
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FIG. 1. Schematic illustration of the atoms with the structure of M 0 S 2 (1H) and Cdl 2 (IT). The 
lattice parameters are a and c. 

II. CALCULATION DETAILS 

Out of the ICSD database, only 16 compounds with MX 2 are semiconductors[T8]. There 
are more layered compounds with more chemical elements. The complexity may hinder their 
practical applications. We selected 14 MX 2 compounds which crystallizes into two different 
hexagonal crystal structures. One is the same as M 0 S 2 and the other is Cdl 2 . They are also 
called 1H and IT structure, respectively. The difference is that the anion hexagonal nets 
are A-A stacked in M 0 S 2 , while they are A-B stacked in Cdl 2 as shown in FigjTJ Thus, the 
first one has a mirror plane of the metals and is lack of inversion, while the second one has 
the inversion center. 

The calculations were performed mainly with pseudo-potential code PWscf|HJj, and cross¬ 
checked by the full-potential local-orbital code (FPLO)[2n| in its current version with the 
default basis settings. All calculations were done within the scalar relativistic approximation. 
The local density approximation (LDA) functional was chosen to be that parameterized by 
Perdew and Wang[21]. The plane-wave kinetic energy cutoff was set to 90 Ry with density 
cut-off of 900 Ry. A shifted 17 x 17 x 2 Monkhorst-Pack mesh was used to perform Brillouin 
zone integration in order to ensure the convergence of the results. Convergence of the total 
energy was set to be better than 10~ 8 Hartree. A vacuum layer with thickness of 30 a.u. 
was used to model the 2D-nature of the compounds. Forces on the atoms are limited within 
0.001 Ry/a b to obtain the ground states and total energy after full geometry relaxation. 
Phonon frequencies and phonon eigenvectors of at the T-point were calculated in a 4 x 4 x 1 
g-grid under the density functional perturbation theory (DFPT). The electron mobilities (//) 
were calculated by summation of the contributions of different scattering sources: 
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where hla is the longitudinal acoustic phonon determined mobility, hop is the longitudinal 
optical phonon and fipz is the piezoelectricity determined ones, respectively. 


III. RESULTS AND DISCUSSIONS 

A. Electron and phonon dispersions 

The electronic structures have been discussed in our previous work|T7]. The most im¬ 
portant parameters related here is the effective masses of electron. As seen in the previous 
works, the compounds with the lT-structure are mostly anisotropic, while the one with 
M 0 S 2 structure are more or less isotropic. The effective electron masses (riT) and the LDA 
bandgap are listed in Table |TJ The m* of M 0 X 2 are about half of that of the free electron 
(m e ). The values of WXo are even lighter which are only about | of m e . With the increase 
of atomic number of X, it slightly increases. 

There is a sizable LDA bandgap in all of these compounds ranging from 0.31 eV in ZrSe 2 
to 1.99 eV in WSe 2 , although it is usually systematically underestimated by LDA. The 
electronic band structure and the atomic projected Density of States (pDOSs) are shown in 
Fig. @ and [3] for the two structures. In the lT-compounds, the conduction band minimum 
is located at the M-point except PtX 2 . It lies between T — M for X = S and Se, and T — K 
for PtTe 2 - There is another local conduction band minimum (Q-valley) only tens of meV 
higher than the global conduction band minimum in the compounds with lH-structure. It 
gives additional electron scattering which will lower the electron mobility. This is also the 
case in PtTe 2 with the lT-structure. The DOS is projected onto different atomic orbital 
orbitals as shown in the figure. The transition metal and the X-atoms both contribute to 
the states of the conduction bands in the 1H compounds. In lT-structure, the situation is 
more complicated. I 11 PtX 2 , where the Pt and X-atoms contribute equally to the states. 
However, in HfX 2 and ZrX 2 , the contributions from the X-atoms are tiny and it is mainly 
due to the metal atoms. They are reversed in SnX 2 . 

The phonon dispersions of these compounds were calculated by the DFPT method. The 
dispersions of the 1H and IT structures are shown in Fig. [4] and [5j The characteristic 
homopolar phonon of the two dimensional monolayer compounds are highlighted in the figure 
with red dashed lines. Within the compounds with the same metal atom, the frequency of 
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TABLE I. The lattice constant, effective mass and the bandgap of the compounds. The effective 
mass in the T — K direction for the M 0 S 2 structure and T — M direction for the Cdl 2 structure is 
calculated. 


mx 2 

a c 

m r-A'(Af) 

m K-M 

EeZ —ph 



(a.u.) (a.u.) 

( m e) 

( m e) 

(eV) 

(eV) 

M 0 S 2 

5.927 2.962 

0.45 

0.45 

3.90 

1.85 

MoSe 2 6.168 3.156 

0.52 

0.52 

3.65 

1.59 

MoTe 2 6.618 3.411 

0.53 

0.57 

0.92 

1.22 

WS 2 

6.047 2.992 

0.24 

0.26 

3.92 

1.99 

WSe 2 

6.166 3.164 

0.33 

0.31 

3.78 

1.71 

SnS 2 

6.879 2.797 

2.11 

0.21 

3.55 

1.42 

SnSe2 

7.165 2.999 

2.91 

0.17 

2.91 

0.65 

HfS 2 

6.731 2.750 

3.30 

0.24 

1.31 

1.05 

HfSe 2 

6.944 2.978 

3.10 

0.18 

1.08 

0.36 

ZrS 2 

6.817 2.771 

1.62 

0.31 

1.52 

1.04 

ZrSe2 

7.007 3.008 

2.03 

0.22 

1.25 

0.31 

PtS 2 

6.670 2.327 

0.26 

0.25 

3.63 

1.69 

PtSe2 

6.978 2.464 

0.21 

0.19 

2.86 

1.25 

PtTe 2 

7.485 2.634 

0.90 

0.77 

1.73 

0.61 


the phonons increases with the atomic number. This is a consequence of the decrease of the 
chemical bond strength, where the forces between the atoms decrease, and the atomic weight 
increases. As the row number of X increases, the atoms become heavier, and the bonds 
become less ionic. Changes of phonon energy are systematic. Taken PtX 2 as an example, the 
width of the acoustical branches decease from 200 cm -1 to 120 cm -1 . The highest optical 
phonon frequency deceases from 380 cm -1 to 180 cm -1 . The gap between the optical and 
the acoustic branches also decreases with the decrease of the atomic weight differences. 
Especially, there is a gap about 100 cm -1 between the optical and acoustic branches in 
PtS 2 - However it diminishes in the other compounds with lT-structure. The gaps also 
diminish in HfSe 2 and ZrX 2 - The larger phonon bandgap is favorite to produce larger 
thermal conductivity due to the reduced acoustical and optical phonon inter-scattering. As 
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FIG. 2. (color online) Electronic bands of 1H structure compounds and the projected and total 
DOS. 


calculated by Gu et a/[22j, the thermal conductivity of lH-type MX 2 are above 50 W/mK 
at room temperature, while that values of the IT-type are much lower. The flat homopolar 
branch are quite obvious as it is a characteristic of these layered structure compounds. The 
linear dispersion of the TA branches becomes parabolic as the atomic number increases in 
PtX 2 due to the rapid decay of the transverse force constants. 

Since the frequency of the optical phonon mode in these compounds is about tens of 
meV, which is comparable with the thermal energy of 30 meV at T= 300 K. It can have 
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FIG. 3. (color online) Electronic bands of IT structure compounds and the projected and total 
DOS. 



















































FIG. 4. (color online) The phonon dispersions of 1H structure compounds. The homopolar mode 
dispersion is in red. 

large populations at room temperatures. Thus both the acoustic and optical phonons may 
have the same order of contribution to the scattering processes. 


B. Acoustic phonon scattering 


The electron mobilities at T = 300 K were approximated by the relaxiation time ap¬ 
proximation. Within the deformation potential approximation [ 231121 ], the electron mobility 
(Takagi model) is approximated by 


Hla 


eh 3 pV s 2 

k B Tm*m d E 2 d _ ph 


( 2 ) 


where kg and h are the Boltzmann constant and the Planck constant, respectively, rn* is 
the effective mass of electron in the electron propagation direction and is the electron 
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FIG. 5. (color online) The phonon dispersions of IT structure compounds. The homopolar mode 
dispersion is in red. 

density of state mass, which is y /m\\m±, where m\\(±_) is the mass parallel (perpendicular) to 
the propagation direction, p is the mass density of the material and V s is the sound velocity 
in the corresponding direction. The validity of this model has already been demonstrated. 
Bruzzone and Fiori used this model to compute the electron mobility of hydrogenated and 
fluorinated graphene as well as h-BCN from first principles[25] and show that graphene with 
a reduced degree of hydrogenation can compete with silicon technology. The sound velocity 
(V s ) is calculated by the superccll method, where frozen phonon mode corresponding to the 
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TABLE II. The sound velocity (V s ), the acoustical(D ac ) and the optical (D op ) deformation poten¬ 
tial. 


mx 2 

v s 

D aC 

Dop 


(km/s) 

(eV) 

(10 8 eV/cm) 

M 0 S 2 

7.93 

3.90,2.4[16],4.5[n] 

1.75,5.8[l6] ,4.1 HI] 

MoSe 2 

6.01 

3.65,3.4[T6] 

1.10, 5.2[T6] 

MoTe 2 

5.04 

0.92 

1.34 

WS 2 

6.67 

3.92,3.2 [16] 

2.34, 3.im 

WSe 2 

5.55 

3.78,3.2[ig 

1.12, 2.3p6] 

SnS 2 

6.18 

3.55 

0.69 

SnSe 2 

4.83 

2.91 

0.38 

HfS 2 

5.86 

1.31 

0.99 

HfSe 2 

4.72 

1.08 

0.62 

ZrS 2 

7.21 

1.52 

1.12 

ZrSe 2 

5.42 

1.25 

0.75 

PtS 2 

6.61 

3.63 

1.06 

PtSe 2 

4.73 

2.86 

0.84 

PtTe 2 

4.89 

1.73 

0.95 


longitudinal phonon with vector q = ^(1,0,0) was simulated. The phonon frequency (a 4 .) 
was obtained, which is related to the sound velocity by 04 = K|q|- The electron phonon 
coupling (E ei-ph) was approximated by the deformation potential D ac . It is related to the 
variation of the electron eigenvalue, which is caused by the volume changes by 

AH 

A E k = D ac —. (3) 

The related values are listed in Table |TT} The values of D ac calculated from fitting the 
relaxion time obtained by Kim|16j and KaasbjergpTJ are also listed in the Table. They 
relative deviation is within 17 %. However, the values of M 0 S 2 obtained by Kaasbjerg is 2.4 
eV, compared with 4.5 eV by Kim and 3.9 eV by us, we can see that this value is method 
sensitive. But the scale in all these calculations is the same. 
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FIG. 6. The different optical phonon modes of 1H and IT structures. 

C. Optical phonon scattering 

The zeroth order optical deformation potential D op is defined as 

A E = D op d (4) 

where A E is the band shift under the proper deformation determined by the proper optical 
phonon mode, d is the atom displacement of the corresponding mode. There are six optical 
modes as shown in Fig. [ 6 j corresponding to six deformation potentials of the conduction 
band minimum. The largest optical deformation potentials are listed in Table [Til The largest 
deformation potential was caused by the A lg (homopolar) mode of the phonon where X 
atoms move in the opposite direction. The only exceptions are PtS 2 and PtSe 2 where the 
contribution are from the E" mode. As reported previous studies|26j. the A lg mode is 
sensitive to electron doping and a red shift is observed on doping, which will increase the 
scattering rate. Thus we can expect a significant decrease of the mobility. However, the 
in-plane mode is less sensitive to doping and the high mobility is expected to be kept. 

The optical phonon limited mobilities are calculated by Equ. (|5]) according to the relax¬ 
ation time approximation: 

_ _ 2 eh 2 pu v _ 

^° P g d m d m*Dl p [N u A 1 + (N u + 1)A 2 ] ’ 

where g d is the valley degeneracy for the final electron states, is the occupation number of 
phonon with angular frequency c o u of mode v which is governed by Bose-Einstein distribution, 
and Ai^ are the onsets of scattering for phonon absorption and emission, both of which are 
set to 1 for intravalley optical phonon scattering. 
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D. Piezoelectric effects 


When the inversion center is missed in the polar semiconductors, optical deformation 
leads to electric polarization which will be a scattering source of carriers. The piezoelectricity 
induced scattering can be modeled in a similar way as the acoustic deformation potential. 
The deformation potential S is replaced by [12] 

r -,2 1 / e n e \2 

2 e r eo 


where eo( r ) is the vacuum (relative) permeability. The piezoelectric coefficient en is defined 
in 2D materials as linear response indicates piezoelectric effect could be treated as the first- 
order coupling between surface polarization and strain tensors (e^). The relation at fixed 
electric held E and temperature T is given by: 

etjk = (7) 

0£jk 

The symmetry required that they are related by 


eni = en 


ei22 = e i2 = -en 


( 8 ) 


e2i2 — e22i — e26 — —en 


With the above relations, the piezoelectric coefficient is calculated by 

-Pi(£n,£22 = 0) — Pi(£u — 0,^22 = 0) = en£n. 


(9) 


In our calculation, the strain, ranging from -0.006 with steps of 0.002 to 0.006, is applied 


to the compounds. The calculated en’s are listed in Table III with comparison with the 
previous work [27]. 

It can be seen that the agreement is good, although our LDA values are systematically 
smaller than the GGA values reported by the work of Duerloo[27]. The en’s of the MoX 2 
compounds are larger than the WX 2 compounds and decrease with the increase of the 
atomic number of X (X = S, Se, Te). The values are in the same order of bulk piezoelectric 
materials which makes significant contributions to the scattering of the electrons. However, 
this scattering can be screened by dielectric effects as shown in experiments. This might be 
the reason where the mobility of the electron can be tuned by dielectric gating. 


13 





TABLE III. The piezoelectric constant (en), dielectric constants (e r ) and the equivalent piezoelec¬ 
tric scattering potential r. pz 


mx 2 

en(10 10 C/m) 

this work ref. [27] 

€ r L 

B p *(eV) 

M 0 S 2 

2.98 

3.06 

4.26 

5.59 

MoSe2 

2.68 

2.80 

4.74 

4.54 

MoTe2 

2.57 

2.98 

5.76 

3.56 

WS 2 

1.72 

2.20 

4.13 

3.32 

WSe 2 

1.51 

1.93 

4.63 

2.61 
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FIG. 7. The total electron mobility and LDA band gap of the 14 compounds with 1H and IT 
structure. 

E. The total mobilities 


The total electron mobilities were calculated by the Matthiessen’s rule as shown in Equ. 
Q, where the contributions from the longitudinal acoustic and optical phonons and the 


piezoelectric scattering were presented in the above sections. The data are listed in Table IV 


It can be seen that the optical phonon largely reduces the electron mobility, which is the main 
source of scattering. We can see that HfSe 2 and ZrSe 2 shown larger acoustic phonon limited 
mobilities, the optical phonon scattering are strong which reduces the electron mobility to 
be less than 100 cm 2 V -1 s -1 . The total mobility and the LDA band gap are shown in Figj7} 
Among the compounds with the 1H structure, WS 2 and WSe 2 have mobilities of 968 and 734 
cm 2 V _1 s _1 , respectively. The LDA band gap is about 1.8 eV, which is only half of that of Si 
(~1750 cm 2 V -1 s -1 ). The well-studied M 0 S 2 has the upper limit of 237 cm 2 V _1 s _1 according 
to our calculation. The experimental value achieved now is about 200 cm 2 V _1 s _1 , which is 
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TABLE IV. The calculated electron mobilities (unit: cm 2 V x s x ) contributed from the different 
scattering sources. 


mx 2 

HLA 

hop 

hPZ 

h 

MoS 2 

1362 

506 

663 

237 

MoSe2 

963 

420 

621 

199 

MoTe 2 10104 

357 

673 

228 

WS 2 

4415 

1553 6148 

968 

WSe 2 

2822 

1191 5921 

734 

HfS 2 

7334 

48 

- 

47 

HfSe 2 

14317 

59 

- 

59 

PtS 2 

4429 

3027 

- 

1795 

PtSe2 

7568 

4746 

- 

2916 

PtTe2 

1467 

235 

- 

202 

SnS 2 

1224 

176 

- 

154 

SnSe 2 

1788 

329 

- 

278 

ZrS 2 

4989 

65 

- 

64 

ZrSe2 

9264 

67 

- 

67 


already quite near this theoretical upper limit. However, in the WS(Se )2 compounds, the 
experimental values are only about tens of cm 2 V _1 s _1 , which is far below our theoretical 
limit. The electron mobilities of Mo compounds are low, which is mainly caused by optical 
phonon scattering. This is in agreement with previous calculation and analyses by Kaasbjerg 
et al. 

In the compounds with lT-structures, the piezoelectric scattering is lacking. Two com¬ 
pounds PtS 2 and PtSe 2 show the upper limits of electron mobilities of 1759 and 2916 
cm 2 V _1 s _1 , respectively. As can be seen from the above sections, the LDA band gaps 
are 1.69 and 1.25 eV. These compounds are thus promising when compared with Si. 

In this work, we have considered only the long wave longitudinal acoustic and optical 
phonon scattering. The initial and final states of the scattered electrons are limited to the 
bottom of the conduction band. That is to say, only the intravalley scattering is included. 
There are other scattering processes, such as interband scattering, and other scattering 
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sources like impurities, electrons and so on. The mobility will be limited by any one of these 
mechanisms. There are uncertainties to precisely determine the contributions from each 
of these mechanism, both experimentally and theoretically, so that we cannot predict the 
mobilities precisely. However, by computing selected scattering sources, the upper bound of 
the mobility can be set. We thus can say, it is hopefully that we can find compounds with 
possible high mobility among the selected ones with larger upper bounds. It is NOT possible 
to find larger mobility in the compounds which are predicted to have lower theoretical upper 
bounds. More sophisticated calculations are needed to aim more precisely within the reduced 
number of compounds. However, these calculations are time consuming and not suitable for 
screening within the large amount of candidates. 

According to our estimation, the mobility of WSe 2 and WS 2 may be larger than M 0 S 2 . 
The mobility order is WS 2 > WSe 2 >MoS 2 >MoSe 2 , which is the same as more elegant 
results by Kim et al. |XS] - Recently, mobilities of WSe 2 and M 0 S 2 are extracted from the 
transfer character curves of field-effect transistors [28]. It is shown that the electron mo¬ 
bility in WSe 2 is about 110 cm 2 V -1 s -1 , while that of MoS 2 is about 25 cm 2 V~ 1 s _1 . These 
experimental results can be a preliminary confirmation of our prediction. 


IV. CONCLUSIONS 

In this work, the electron mobility of 14 MX 2 type two dimensional semiconductors were 
calculated where only elastic scattering from long wave acoustic and optical phonons were 
taken into account by the deformation potential approximation. The piezoelectric scattering 
was included in compounds lack of inversion center and treated as the acoustic phonon-like 
deformation potential. We found that the total electron mobility in WS 2 can reach 968 
cm 2 V -1 s _1 , which is the highest value among the compounds with the lH-structures. The 
value of 2916 cm 2 V~ 1 s~ 1 can be reached in PtSe 2 with the lT-structure, where the inversion 
center is presented and the piezoelectric scattering is prohibited. WS 2 is a direct band 
gap semiconductor with band gap of 1.99 eV at the K-point, while PtSe 2 is an indirect 
band gap of 1.25 eV from the T-point to the middle of P-M line. Concerning with these 
two requirements, theses two compounds are promising for two-dimensional semiconductors 
used for future logical devices. 
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